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1.  Background 

The  assignment  came  about  with  the  main  aim  of  developing  new  thrusts  for  the  Air 
Force  Research  Laboratory,  Sensors  Directorate,  Radar  Signal  Processing  Branch 
(SNRT)  with  strong  ties  to  security  applications  such  as  perimeter  surveillance,  concealed 
weapon  detection  (CWD)  and  related  problems.  My  own  previous  work  in  image 
processing  for  CWD  and  other  security-related  imaging  with  visual,  x-ray,  infrared  and 
millimeter  wave  imagery  was  seen  as  a  jumping-off  point  for  identifying  new  areas.  The 
IP  A  agreement  listed  the  following  as  anticipated  roles: 

1.  To  define  key  research  problems  in  homeland  security  that  are  amenable  to 
solution  in  a  scientific  framework. 

2.  To  initiate  a  research  program  to  investigate  all  or  a  subset  of  the  problems 
identified. 

3.  To  interface  with  the  Department  of  Homeland  Security  and  other  agencies  for 
transitioning  of  research  product  for  actual  use. 

4.  To  ensure  that  the  program  becomes  self-sustaining  with  potential  to  contribute 
over  the  next  decade. 

5.  Items  1  and  2  have  largely  been  fulfilled  and  efforts  are  under  way  in  regards  to 
items  3  and  4. 
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2.  Identification  of  Opportunity 

A  substantial  part  of  the  early  effort  was  devoted  to  identifying  the  most  appropriate  new 
thrust.  Several  different  possibilities  were  considered  including  sensor  networks,  new 
evolutionary  programming  based  optimization  techniques  and  signal  processing  for 
concealed  object  detection  and  classification.  A  detailed  evaluation  of  the  state  of  the  art 
of  these  different  fields  was  performed.  New  innovations  in  the  first  two  fields  are  much 
more  tied  to  generic  communication  and  architectural  problems  than  to  a  direct 
addressing  of  the  problem  of  detecting  and  classifying  threats.  On  the  other  hand,  the 
third  area  addresses  the  security  and  surveillance  problem  directly  and  is  well  in  line  with 
SNRT’s  technical  interests. 

We  have  previously  addressed  the  problem  of  CWD  in  the  context  of  detecting  weapons 
concealed  under  clothing.  Infrared  and  millimeter  wave  images  were  considered  for  the 
purpose.  Infrared  detection  depends  on  the  concealed  weapon’s  temperature  being 
significantly  different  from  that  of  its  surroundings.  Millimeter  wave  (MMW)  detection 
relies  on  reflection  from  the  object  of  interest.  Passive  detection  is  possible  in  many 
instances  since  MMW  is  present  naturally  and  can  penetrate  various  types  of  clothing. 
Active  devices  can  also  be  used  to  provide  MMW  illumination.  Examples  of  our  prior 
work  in  enhancing  such  images  are  found  in  [1].  It  was  found  that  Infrared  (IR)  images 
are  less  reliable  since  one  can  easily  defeat  or  mask  the  object’s  response.  MMW  detector 
technology  has  only  now  advanced  to  the  point  where  real  time  images  are  possible. 
Imaging  in  the  Terahertz  (THz)  wavelengths  was  an  avenue  that  suggested  itself  as 
worthy  of  exploration.  Like  MMW  it  penetrates  clothing.  However,  with  wavelengths 
ranging  from  1  micron  to  1  mm,  it  provides  superior  resolution.  Furthermore,  it  offers 
spectroscopic  capability  that  enables  object  identification  not  only  by  shape  attributes  but 
also  by  material  composition  [2],  This  opens  up  applications  such  as  standoff  detection  of 
explosives  [3].  The  range  of  300  GHz  to  3  THz  has  been  explored  by  many  researchers. 
Standoff  ranges  of  up  to  50  meters  can  be  achieved  at  these  frequencies.  The  main 
limitation  is  water  vapor  absorption.  There  is  huge  interest  in  the  potential  of  THz 
technology.  A  considerable  amount  of  work  has  been  done  in  developing  THz  generators 
and  detectors.  However,  there  is  much  work  that  needs  to  be  done  in  the  area  of  signal 
and  image  processing  algorithms  for  making  THz  technology  useful  in  applications  such 
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as  CWD  and  explosives  detection.  The  next  section  will  provide  an  overview  of  THz 
imaging. 
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3.  THz  Technologies 

Although  a  precise  definition  of  the  Terahertz  region  does  not  exist,  a  commonly 
accepted  range  is  the  300  GHz  to  3  THz  region  [4].  This  corresponds  to  wavelengths  of 
100  pm  to  1  mm.  This  has  been  a  relatively  unexplored  region  of  the  electromagnetic 
spectrum  and  has  been  referred  to  as  the  “Terahertz  gap.”  The  gap  has  had  mostly  to  do 
with  the  lack  of  source  and  detector  technologies  until  recently.  For  Electro-Magnetic 
(EM)  radiation  with  wavelengths  longer  than  the  THz,  antenna  and  waveguide 
technologies  have  been  available  for  several  decades  whereas  for  shorter  than  THz 
wavelengths,  optica, 1  and  optoelectronic  technologies  have  been  used.  The  two 
technologies  have  to  be  brought  together  to  generate  and  study  the  THz  region,  which  is 
also  referred  to  as  the  far  infrared. 

Some  of  the  earliest  works  that  invoke  the  tenn  Terahertz  are  in  acoustics  in  the  context 
of  photon  generation  [5,  6].  The  occurrence  of  THz  bandwidths  in  the  context  of  high 
pressure  CO2  laser  operation  was  predicted  in  [7],  Explicit  reference  to  a  device  as  a 
potential  detector  of  THz  radiation  occurs  in  [8].  Synthesis  of  terahertz  radiation  has 
been  made  possible  by  the  development  of  three  distinct  approaches  [9,  10]: 

1 .  Direct  emission  from  far-infrared  incoherent  sources, 

2.  Up-conversion  from  microwave  frequencies  to  terahertz  frequencies. 

3.  Down-conversion  of  coherent  lasers  operating  in  the  near  infrared  to  the  terahertz 
band. 

The  optically  driven  sources  exhibit  more  power  at  terahertz  frequencies  than  direct 
emission  or  up-conversion.  A  typical  arrangement  for  pulsed  THz  generation  is  shown  in 
the  context  of  time  domain  spectroscopy  as  shown  in  Figure  1.  A  light  pulse  from  the 
femtosecond  laser  source  strikes  a  device  (shown  separately  in  Figure  2)  which 
essentially  consists  of  electrodes  embedded  in  a  Gallium  Arsenide  (GaAs)  substrate.  A 
Direct  Current  (DC)  bias  is  maintained  across  the  electrodes.  The  current  discharge  due 
to  ionization  of  the  dielectric  from  the  laser  pulse  results  in  the  generation  of  a  free  space 
THz  pulse.  The  same  device  works  as  a  detector  in  reverse.  Instead  of  maintaining  a  DC 
bias,  wires  are  led  away  from  the  electrodes.  When  a  THz  pulse  strikes  the  substrate  in 
synchronization  with  a  light  pulse,  a  current  pulse  is  generated  in  the  wires. 
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Figure  1.  THz  time-domain  spectroscopy  using  pulsed  THz  generation.  After  Eric  Mueller  [10]. 
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Figure  2.  THz  generation  device.  After  Eric  Mueller  [10]. 


The  attributes  of  various  methods  of  THz  generation  are  summarized  in  Table  1  [10]. 
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Table  1.  Properties  of  THz  generation  methods 


Generation  Method 

Optically 

pumped  THz 

lasers 

Time- 

domain 

spectroscopy 

Direct 

multiplied 

sources 

Frequency 

mixing 

Attribute 

Average 

power 

>  100  mW 

~  1  pW 

mW  -  pW 

<  1  pW 

Frequency 

range 

0.3  to  10  THz 

0.1  to  2  THz 

01.  To  1  THz 

0.3  to  10  THz 

Tunability 

Discrete  steps 

N/A 

10  to  15%  of 

center 

Continuous 

Mode 

CW  /  Pulsed 

Pulsed 

CW 

CW 

3.1  Time-Domain  Spectroscopy 

The  arrangement  in  Figure  1  can  be  used  to  perform  time-domain  spectroscopy  (TDS). 
The  optical  delay  shown  in  the  figure  is  used  to  adjust  the  time  instant  of  sampling  THz 
radiation  from  the  detector.  The  material  sample,  whose  THz  spectral  response  is  to  be 
obtained,  is  positioned  in  the  path  of  the  generated  THz  pulse  and  the  output  is  sampled 
at  the  detector  at  different  instants.  If  the  THz  pulses  are  generated  in  periodic  fashion, 
the  optical  delay  can  be  adjusted  to  gate  the  detector  at  different  instants  in  the  period. 
The  samples  at  the  various  instants  can  be  used  to  construct  a  time-domain  profile  of  the 
response  of  the  material.  A  Fourier  transformation  of  the  profile  provides  the  spectral 
response  of  the  material  in  the  THz  regime.  An  image  can  be  generated  by  optical 
scanning  of  the  input  beam  or  by  moving  the  sample  around.  THz  TDS  is  by  no  means 
confined  to  transmission  data.  It  is  possible  to  get  reflective  data  from  samples  by  slight 
alteration  of  the  detector  path  in  Figure  1. 

3.2  Terahertz  Tomography 

The  fact  that  full  field  measurements,  that  is,  phase  and  magnitude  measurements  are 
available  with  THz  makes  reflection  tomography  possible  [11].  In  principle,  range 
calculations  are  made  by  comparing  the  phase  of  THz  reflections  with  the  phase  of  the 
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transmitted  THz  signal.  The  range  map  as  a  function  of  position  allows  for  reconstruction 
of  the  surface  of  the  object  in  three  dimensional  (3D).  The  range  measurements  can  be 
made  with  either  pulsed  or  continuous  wave  THz  waveforms.  A  version  of  tomographic 
reconstruction  resembling  seismic  migration  is  proposed  in  [12].  As  shown  in  Figure  3, 
pulses  sent  from  transmitter  T  are  received  after  reflection  at  layer  boundaries  by  multiple 
receiver  elements  R.  Delay  and  sum  operations  are  performed  to  reconstruct  reflection 
intensity  at  various  points. 


RRRRRRRRR  T  RRRRRRRRR 
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4.  Perspective  and  Summary  of  Accomplishments 

Figure  4  summarizes  our  perspective  of  the  components  required  for  successful 
development  of  THz  imaging.  It  involves  the  confluence  of  three  entities:  signal 
processing,  photonics,  and  electronics. 


Figure  4.  Components  of  THz  imaging 

Clearly  a  lot  of  effort  has  been  invested  in  the  photonics  and  electronics  areas  through  the 
development  of  generation  and  detection  technologies.  The  area  that  is  least  developed  is 
signal  processing  of  TFIz  data.  In  this  regard,  the  key  accomplishments  of  the  assignment 
have  been  the  following: 

1.  Key  signal  processing  issues  that  need  to  be  addressed  for  significant 
improvement  in  THz  imaging  have  been  identified  and  potential  lines  of 
investigation  have  been  proposed. 

2.  Preliminary  mathematical  formulation  and  simulation  results  have  been  obtained 
for  signal  processing  problems  in  THz  image  reconstruction. 

Focal  plane  detector  arrays  are  not  yet  available  in  the  usual  operating  range  of  0.3  to  5 
THz.  Antenna  array  based  imaging  using  signal  processing  techniques  is  therefore  an 
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attractive  alternative.  We  have  developed  mathematical  formulation  and  preliminary 
solutions  for  correlative  interferometric  image  reconstruction  problems  with  antenna 
arrays.  Computer  simulations  have  also  been  performed  for  tomographic  reconstruction 
of  complex-valued  projections  using  MATLAB  phantoms.  The  next  section  details  the 
investigation  into  correlative  interferometry. 
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5.  Imaging  through  Correlative  Interferometry 

The  construction  of  THz  images  of  extended  objects  using  correlative  interferometry  has 
been  proposed  by  Federici  and  colleagues  [13].  The  key  advantages  of  this  approach  are: 

1 .  It  works  even  with  incoherent  radiation. 

2.  Synthetic  aperture  techniques  can  be  used  to  reduce  the  number  of  actual 
detectors  used. 

Researchers  in  this  field  suggest  that  improvement  to  the  synthesized  data  can  be  done 
with  techniques  such  as  the  CLEAN  algorithm  and  the  Maximum  Entropy  Method 
(MEM)  that  have  been  widely  used  in  radio  astronomy.  Our  view  is  that  such  techniques 
are  oriented  towards  point  sources  and  scatterings.  Furthermore,  the  CLEAN  algorithm 
requires  at  least  the  first  identified  peak  to  be  a  true  fully  resolved  isolated  component. 
We  propose  developing  constrained  optimization  methods  that  offer  super-resolution 
potential. 

5.1  Correlative  interferometry  in  the  far-field 

Suppose  /(0)cos[27r/t  +  cp(0)]  is  the  far-field  electric  field  amplitude  at  x  -  0  of  the 
EM  wave  with  direction  of  arrival  0  (see  Figure  5). 


Extended  object 
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Then,  the  contribution  of  the  field  covering  an  elemental  angle  c/0  in  the  neighborhood 
of  this  angle  to  the  field  at  x  is  equal  to 


/(9)cos  2nf^t~  — sin0j  +  (p(9)  d 9.  (1) 

The  total  field  F(x)at  x  is  equal  to  an  integration  of  the  above  over  all  angles  and  is 
given  by 


F{x)=  J  /(0)cos  2t if  t — sin9  +(p(0)  c/9 


=  [  /(0)cos  27r/zl-^^-sin9  +  (p(9)  c/0 

-in  L  ^ 

The  cross-correlation  between  the  fields  at  two  distinct  points  x  and  y  is  given 


R(x,y)  =  F{F(x)F(y)}  = 


r\ 

e\  |  |  /(O)eos  2  nft — — sin9  +  (p(0)  /(p)cos  27t//--^^-sinp-i-(p(|u)  c/0  c/pi 


2  ny  . 


Equation  (3)  can  be  written  as 


where 


£{F(x)F(y)}=i±i 


Il=E< j  J  J  /(9)/(ju)cos  ^^sin9-^-sinp  +  cp(|j,)-(p(0)  c/9c/p 


I2=e\  J  J  /(9)/(ju)cos  4Tc//--^^-sin9--^-sinp  +  (p(9)  +  (p(iu)  c/9c/pl. 
{-tn-in  l  X  X  J  J 

If  we  have  incoherent  radiation  and  assume  that  cp(0)  is  uniformly  distributed  between 
-7i  and  re  with  cp(0)  and  cp(ju.)  independent  for0  ^  p  ,  we  have 


71/  Z 

R{  x,y)  =  /,  I  v.  i  )  =  J 


sin  9  c/0  . 
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Equation  (5)  shows  that  R(x,y)  is  a  symmetric  function  solely  of  the  displacement 

(x-y) .  In  fact,  it  is  a  true  autocorrelation  function.  The  above  expression  essentially 

relates  the  source  intensity  f2  (9)  to  the  autocorrelation.  A  Fourier  transform 

relationship  can  be  developed  as  follows. 

Let 


>Jx-yX- 

(6) 

Then  the  autocorrelation  relationship  becomes 

n/2  1-2  / qS 

f?(s)  =  J  - — —  cos [2ns sin 6) c/6 

-n/2  ^ 

(7) 

Put 

t  =  sin  0 

(8) 

We  then  have 

00 

^(5)=  J  g(t)cos(27r5t)fi?t 

(9) 

-oo 


Where 

/2(sin_1t) 

g(0=<  (io) 

0  otherwise. 

Equation  (9)  indicates  that  the  autocorrelation  A* (.s') and  the  derived  field  intensity 

function  g(t)  form  a  Fourier  transform  pair.  There  is  thus  the  possibility  that,  from  an 

estimate  of  the  autocorrelation,  the  extended  object  can  be  synthesized. 

Simulations  were  performed  with  a  discretized  version  of  the  problem  by  implementing 
the  relationship  in(7).  Figure  6  shows  an  input  object  (essentially  the  squared  magnitude 
of  field  strength)  generated  as  a  function  of  0  with  243  samples  between  -n  /  2  and  n/2. 


12 


Input  object 


Figure  6.  Input  object. 


Figure  7  shows  243  samples  of  the  autocorrelation  obtained  from  (7)  for  an  array  spacing 
of  0.2  wavelengths.  The  corresponding  243x243  matrix  is  singular  with  a  rank  of  64. 
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Autocorrelation 


Figure  8  shows  the  reconstruction  obtained  from  an  unconstrained  minimization  as 

II  m2 

f  =  argmin  Mx-r  (11) 

x  H 

where  A  is  the  matrix  transforming  the  object  to  the  autocorrelation  through  the 
discretization  of  (7)  and  r  is  the  recorded  autocorrelation.  In  practice,  the  autocorrelation 
is  not  available  and  has  to  be  estimated  from  data  samples.  The  unconstrained 
reconstruction  in  Figure  8  is  grossly  different  from  the  true  object. 
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Unconstrained  Reconstruction 


Figure  8.  Result  of  unconstrained  optimization 


Imposing  a  positivity  constraint  on  the  solution  yields  the  reconstruction  shown  in  Figure 
9. 
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Intensity 


Reconstruction  with  positivity  constraint 

1.5 


1 


0.5 
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A  noisy  reconstruction  of  the  taller  component  of  the  object  seems  to  result  from  this 
constraint.  Limiting  the  reconstruction  boundary  to  the  true  limits  results  in  the 
reconstruction  of  shown  in  Figure  10. 


e 

Figure  9.  Result  of  applying  positivity  constraint  only. 
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Intensity 


Reconstruction  with  boundary  constraint 


Figure  10.  Result  of  applying  boundary  constraint  only 

The  positivity  constraint  was  not  imposed  and  it  is  seen  to  be  violated.  Figure  1 1  shows  a 
reconstruction  with  both  the  positivity  and  boundary  constraints  applied. 


17 


Intensity 


Reconstruction  with  positivity  and  boundary  constraints 

1.2 
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0.8 

0.6 

0.4 
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The  result  is  obviously  a  much  closer  reconstruction  to  the  original.  Figures  12  and  13 
show  reconstructions  at  two  different  SNRs. 


1.5  -1  -0.5  0  0.5  1  1.5 

0 

Figure  11.  Both  constraints  applied 
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Intensity 


Reconstruction  at  3  dB  SNR 


Figure  12.  Reconstruction  in  noise.  SNR=3  dB.  Both  constraints  applied. 
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Reconstruction  at  0.6  dB  SNR 


Figure  13.  Reconstruction  in  noise.  SNR=0.6  dB.  Both  constraints  applied. 


5.2  Near-Field  Imaging 

Terahertz  imaging  with  antenna  arrays  has  been  proposed  by  Federici  and  colleagues 
[14].  They  indicate  that  high  resolution  imaging  necessarily  involves  working  in  the  near¬ 
field.  This  opens  the  substantially  unexplored  subject  of  array  signal  processing  for  near¬ 
field  data.  This  document  contains  preliminary  results  on  the  modeling  and  simulation 
using  one-dimensional  fields. 

The  field  at  a  distance  r  from  a  point  source  field  of  strength  E  is 

£ 

—  expy'(o0(t-r/c)  (12) 

r 


This  may  be  rewritten  as 


—  exp  j  (co0t-27tr/A,) 

where  X  is  the  wavelength.  Thus  the  field  phasor  changes  by  an  amount 


(13) 
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(14) 


—  exp  (-j  2  nr /A,) 

at  a  distance  r  from  the  source.  Consider  a  line  source  and  another  line  parallel  to  it  at 
distance  D  wavelengths  from  it  as  shown  in  Figure  14. 


Source  Sensing 

line 

Figure  14.  Source-sensing  geometry 


If  E  (x)  is  the  source  field  phasor  at  position  xk,  the  contribution  from  an  elemental  unit 
of  length  d[xX)  in  its  neighborhood  to  the  field  at  position  yX  on  the  sensing  line  is  given 
(using  Eq.  (14))  by 


E(x) 


^exp \-j2nJD2  +(y-x)  \dx 
lD2+(y-x)z  1  '  ' 

Thus,  the  sensed  field  S'(y)  at  yX  is  expressed  as  the  convolution 

S(y)  =  E(y)*h(y) 


(15) 


(16) 


where  *  denotes  the  convolution  operation,  and 
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exp 


/z(x)  =  A- 


r~ 

f  v  \ 

2 

-j2n  l+\ 

[  X  ' 

V  1 

\Dj 

1  + 


-  with 


\Dj 


A  = 


cxp(-/27iD) 


D 


(17) 


The  term  X  cancels  out  as  a  common  factor.  Equation  (16)  approximates  to  a  Fourier 
transfonn  relationship  in  the  far  field  [15].  However,  our  interest  is  in  working  with  this 
equation  as  is  since  we  are  considering  the  near-field.  The  function  h  (x)  may  be  regarded 


as  an  impulse  response  of  a  linear  time-invariant  system  or  as  a  Green’s  function  for  the 
particular  geometry. 

Suppose  we  have  N  observations  at  sv...,sN  at  positions  y  =  Xxl,..., XxN  respectively  on 


the  sensing  line.  We  address  the  problem  of  reconstructing  the  source  field  E (x)  from 


these  observations.  From  (16),  we  have 

sn=E(xyh(x)\x=x* 

=  ^E(\)h(xn-\)dv,  n  =  l,...,N. 


(18) 


The  problem  is  one  of  reconstructing  E  (x)  from  the  known  function  h  (x)  and  the  values 
of  sn  andx„.  Obviously,  it  is  generally  impossible  to  reconstruct  a  continuous  function 
from  a  finite  number  of  equations.  One  may  find  a  least  squares  solution  by  finding 
E (x)  that  minimizes 


n 

This  will  mean  finding  a  solution  to  the  equation 

]E(v)i^h(n-x)ti  (n-\i)dv  =^jsnh*  (n-\i)  forallp 


(19) 


(20) 


Variations  of  this  equation  appear  in  various  sampling  contexts  [16].  An  expression  can 
be  obtained  in  tenns  of  the  Fourier  transfonn  off? (x) .  However,  it  is  unwieldy  and  does 
not  point  to  how  one  might  solve  (20).  Furthermore,  it  should  be  clear  that  there  are  many 
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potential  solutions.  Further  constraints  need  to  be  imposed  to  generate  a  unique  solution. 
We  develop  the  approach  through  simulation  examples  of  a  discretization  of  the  problem. 
Correlative  interferometry  can  be  developed  for  the  near  field  as  well.  Suppose  we  now 
model  the  source  field  at  point  xA  as 

E  [x,t)  —  E  (x)exp  jui0t  —  |  E  (x)lexp  j  {u0t  +  <bv)  (21) 


where  <hvis  a  random  variable  uniformly  distributed  between  -7r  and  n .  We  assume 
further  that  this  phase  tenn  is  spatially  uncorrelated,  that  is, 

£{*,*,}  =  0  (22) 

for  x^y .  Here  E{  }is  the  expectation  operator.  Equation  (22)  corresponds  to  situation 

where  we  have  a  non-coherent  source  or  illumination. 

The  spatial  correlation  of  the  sensed  field  is  given  by 

\E{l^)\f{x~X)f(y-n)  (23) 

x  E  {exp  j  (<bA  —  )|  dX  d/i 

where 


fix) 


exp 


~JK 


yD, 


1  + 


x 

D 


(24) 


It  follows  from  (22)  that 

£  {exp  7  (<bA— $„)}  =  ()  for  A  ^ /x.  (25) 

Also,  the  right-hand  side  of  (23)  is  not  a  function  of  t.  Therefore,  we  have 
f?(x,y)  =  (x,t)S*  (>>,*)} 

=  -F  n£(A)|7(Jp_A)/-(^-A)rfA  (26) 

D~l 

=  j\E[Xf  g[x-\,y-\)d\ 
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where  R{x,y)  denotes  the  spatial  correlation  of  the  sensed  field.  Equation  (26)  shows  that 
it  is  a  linear  transfonnation  of  the  source  intensity  through  the  spatially  varying  kernel 

g[x,y)=  Di  (2y) 

If  R(x,y)  is  available,  then  reconstructing  the  source  intensity |£(A)|"  from  (26)  requires 

generating  a  solution  to  an  integral  equation.  Since  real-world  measurements  are  from  an 
array  that  samples  on  the  sensing  line,  we  turn  our  attention  to  the  problem  of 
reconstructing  with  such  samples.  Suppose  we  are  given  N  spatial  observations 
Sj  ^(t)  at  positions  y  =  Xyl,...,XyN  respectively  on  the  sensing  line.  From  (26),  we 

have 

R  (y  m  >y»)  =  f  \E  (A)|  0 y«  -  A,  y«  ~  A) dX  (28) 

L 

Solving  for  the  intensity  in  (28),  given  the  spatial  correlation  values  between  the  sensors, 
amounts  to  reconstructing  a  continuous  function  from  a  finite  number  of  discrete  values, 
clearly  an  ill-posed  problem.  Further  constraints  are  required  to  generate  a  unique 
solution. 

We  can  discretize  the  problem  further  by  approximating  (28)  in  terms  of  samples  of  the 
source  image  as 

M  2 

R  (. ym ,  y„ )  =  \E  ixk  )|  s  ( ym  -  xk » y«  -  xk )  (29) 

k= 1 

where  M  samples  of  the  source  image  have  been  considered  at  xx,x2,...,xM  respectively.  If 
M  —  N2 ,  one  could  try  solving  (29)  with  as  many  equations  as  unknowns.  If  M  <N2 ,  that 
is,  with  more  equations  than  unknowns,  we  can  use  a  pseudo-inverse  solution.  However, 
the  following  issues  need  to  be  addressed: 

1 .  Constraints  may  have  to  be  imposed  on  the  solution.  At  a  minimum,  a  positivity 
constraint  may  be  needed  since  we  are  reconstructing  the  source  intensity. 

2.  In  practice,  we  have  access  to  the  field  measurements  and  not  to  the  spatial 
correlation. 


24 


We  therefore  propose  a  constrained  minimization  approach  using  estimates  of  the 
correlation.  Suppose  R(ym,y„)i  m,n  =  are  the  estimates  of  the  spatial  correlation. 
We  solve 

argminC(4)  =  ]r]r 

m= 1  n= 1 

subject  to: 

1 .  Positivity  constraint  on  Ik . 

2.  Boundary  constraint  on  Ik . 

Here  the  I  k  s  are  estimates  of  the  source  intensity  at  xt,x2,...,xM  . 

Simulations  were  performed  assuming  availability  of  true  correlation  values.  Figure  15 
shows  an  image  whose  cross-sectional  profile  is  the  one-dimensional  object  of  interest. 


M 


R  ( ym  >  y„ )  -  J2 1  kS  ( ym  -  x * » y«  -  xk ) 


k= 1 


(30) 
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Figure  15.  Input  object 


The  horizontal  axis  consists  of  4 1  samples  taken  over  a  distance  of  400  wavelength  units 
although  the  object  itself  is  much  smaller.  The  41x41  matrix  transforming  the  input  field 
to  the  sensed  field  is  singular.  Figure  16  shows  an  unconstrained  reconstruction  while 
Figure  17  shows  a  reconstruction  obtained  with  a  positivity  constraint  and  a  support 
constraint  which  restricts  the  reconstruction  to  the  support  region  of  the  input  object.  The 
use  of  such  constraints  is  characteristic  of  various  super-resolution  techniques.  Needless 
to  say,  the  support  constraint  results  in  a  reconstruction  closer  to  the  true  field.  Flowever, 
the  difference  clearly  shows  that  further  constraints  may  be  needed  to  generate  a  unique 
solution. 
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Figure  16.  Unconstrained  reconstruction 
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Figure  17.  Reconstruction  with  positivity  and  boundary  constraints. 
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6.  Conclusions  and  Future  Work 

The  assignment  has  identified  Terahertz  imaging  for  standoff  detection  of  concealed 
objects  as  worthy  of  pursuit  by  AFRL/SNRT.  The  future  success  of  THz  imaging  rests 
not  only  on  progress  in  optoelectronic  and  electromagnetic  generators  and  detectors  but 
also  on  the  development  of  signal  processing  algorithms.  Our  investigation  has  shown 
that  interesting  and  important  problems  for  THz  image  reconstruction  from  multisensor 
arrays  can  be  formulated.  Preliminary  solutions  have  also  been  provided.  To  this  date, 
two  publications  have  resulted  from  this  work  Further  investigation  is  necessary  to  test 
the  performance  of  the  algorithms  and  refine  them  for  speed  and  robustness  with  respect 
to  impairments  such  as  noise,  sensor  positioning  errors  etc. 

Some  investigation  was  also  conducted  into  tomographic  reconstruction  of  complex 
fields  in  collaboration  with  2nd  Lt.  Joshua  Markow,  AFRL/SNRT.  The  idea  here  is  to  take 
advantage  of  the  fact  that,  unlike  x-rays  which  offer  only  magnitude  information,  THz 
offers  phase  information,  as  well.  While  the  magnitude  contains  attenuation  information, 
the  phase  contains  infonnation  about  the  refractive  index.  The  investigation  has  revealed 
that  phase  reconstruction  is  very  sensitive  to  noise  and  aliasing.  Future  work  can 
investigate  methods  of  mitigating  the  amount  of  phase  errors. 

A  range  of  THz  frequencies  that  is  yet  to  be  explored  is  in  the  neighborhood  of  30  THz 
(wavelengths  in  the  neighborhood  of  10.6  micron  or  the  carbon  dioxide  wavelength).  As 
pointed  out  by  Domenico  Luisi,  AFRL/SNRD,  there  is  very  little  water  vapor  absorption 
at  these  wavelengths  and  the  radiation  can  travel  much  larger  distances  than  the 
commonly  used  range  that  falls  under  3  THz.  There  is  thus  the  possibility  that  detection 
can  be  accomplished  at  greater  standoff  distances.  However,  substantial  spectroscopic 
work  is  needed  in  this  area. 

Investing  in  the  development  of  new  signal  processing  algorithms  offers  the  greatest  near 
term  (3  to  5  year  time  frame)  return  is  in  interferometric  tomographic  reconstruction  from 
antenna  arrays.  The  key  mathematical  problems  here  are  that  of: 

1.  Relating  field  or  intensity  at  source  to  antenna  array  measurements  at  a  given 
wavelength. 

2.  Estimating  source  field  or  intensity  from  above  relationship  in  the  presence  of 
impairments  such  as  finite  apertures,  illumination  issues  and  noise. 
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The  interferometric  approach  relies  on  exploiting  correlation  and  spatial  diversity  in  the 
antennas  located  at  different  geographic  locations.  Several  different  scenarios  for 
tomographic  reconstruction  in  the  order  of  increasing  complexity  are  outlined  below: 

1 .  Far  field  passive  reconstruction. 

2.  Near  field  passive  image  reconstruction. 

3.  Multistatic  image  construction  with  active  illumination. 

All  of  the  scenarios  above  are  worthy  of  investigation.  The  results  of  scenarios  1  and  2 
are  needed  to  carry  out  investigation  of  scenario  3.  Passive  reconstruction  assumes 
continuous  and  uniform  illumination.  The  reconstruction  algorithms  developed  for 
passive  scenarios  can  be  embedded  as  components  within  the  active  multistatic  image 
reconstruction  algorithms. 
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structure,  characteristics,  arithmetic  and  feasibility  are  discussed.  System  performance  and 
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industry  application  perspective  are  analyzed,  specially  compared  with  X-ray  process 
tomography  system. 

5.  Gregory,  I.S.;  Tribe,  W.R.;  Cole,  B.E.;  Baker,  C.;  Evans,  M.J.;  Bradley,  I.V.;  Linfield,  E.H.; 
Davies,  A.G.;  Missous,  M.;  “Phase  sensitive  continuous-wave  THz  imaging  using  diode 
lasers”  Electronics  Letters  ,  Volume:  40  ,  Issue:  2  ,  Pages:  143  -  145,  22  Jan.  2004. 

Presents  the  preliminary  results  of  a  continuous-wave  all-optoelectronic  terahertz  imaging 
system  based  on  diode  lasers.  The  coherent  detection  scheme  is  phase  sensitive  and  operates 
at  room  temperature.  The  experiment  results  in  low-cost,  compact  system  having  image 
capture  rates  comparable  with  those  from  state-of  the-art  pulsed  THz  systems.  This  is  the  first 
demonstration  of  such  a  compact,  robust  and  low-cost  system. 

6.  Ferguson,  B.;  Wang,  S.;  Xi,  J.;  Gray,  D.;  Abbott,  D.;  Zhang,  X.-C.;  “Linearized  inverse 
scattering  for  three  dimensional  terahertz  imaging”  Conference  on  Lasers  and  Electro- 
Optics,  2003.  CLEO  '03.,  1-6  June  2003. 

Uses  two  dimensional  electro-optic  sampling  that  allows  a  target's  scattered  terahertz  field  to 
be  measured.  By  obtaining  multiple  projection  images,  reconstruction  of  three  dimensional 
targets  with  millimeter  resolution  was  performed.  The  reconstruction  was  performed  using  a 
frequency  range  of  0.2  THz,  which  provided  the  maximum  SNR  for  the  antenna  source  used. 
This  reconstruction  algorithm  is  only  applicable  to  weakly  scattering  objects. 

7.  Buma,  T.;  Norris,  T.B.;  “Time  reversal  3-D  imaging  using  single-cycle  terahertz  pulses” 
Conference  on  Lasers  and  Electro-Optics,  2003.  CLEO  '03.,  1-6  June  2003 

Explains  demonstration  of  time  reversal  3-D  imaging  using  single-cycle  terahertz  pulses.  The 
photoconductive  transmitter  and  receiver  are  arranged  in  a  reflection  imaging  configuration. 
The  synthesized  ring  annular  detection  array  can  distinguish  spatially  separated  image  planes. 

8.  Zhong,  H.;  Xi,  J.;  Wang,  S.;  Zhang,  X.-C.;  “T-ray  tomography  with  a  Fresnel  lens” 
Conference  on  Lasers  and  Electro-Optics,  2003.  CLEO  '03.,  1-6  June  2003 

This  article  details  in  brief  a  demonstration  of  tomographic  imaging  with  Fresnel  lens,  since 
Fresnel  lens  is  lighter  and  compact  than  conventional  THz  optics.  The  article  also  reports 
about  the  possibilities  of  long  range  terahertz  imaging  with  the  same  lens. 

9.  Wang,  S.;  Ferguson,  B.;  Zhong,  H.;  Zhang,  X.-C.;  “Three-dimensional  terahertz 
holography”  Conference  on  Lasers  and  Electro-Optics,  2003.  CLEO  '03.,  1-6  June  2003  . 

This  article  demonstrates  three-dimensional  terahertz  (THz)  holography  with  the  use  of  the 
windowed  Fourier  transformation.  This  paper  explains  that  the  coherent  nature  of  THz  time 
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domain  spectroscopy  allows  measurement  of  the  phase  and  amplitude  of  the  THz  wave  for 
holographic  imaging.  Results  of  the  demonstration  have  not  been  explained  or  shown  in  this 
article. 

10.  Evans-Pughe,  C.;  “Sharper  sight  for  surgeons”  IEE  Review  ,  Volume:  49  ,  Issue:  9  ,  Sept. 
2003  Pages:40  -  43 

This  article  states  that  terahertz  radiation  could  become  a  new  weapon  in  the  fight  against 
cancer  due  to  its  capacity  to  identify  minute  variations  in  soft  tissue  undetectable  by  other 
imaging  methods.  The  radiation  excites  interactions  between  adjacent  molecules  and  it's  the 
secondary  radiation  produced  by  these  molecular  vibrations  that  provides  the  basis  for  cancer 
detection.  The  resulting  spectral  'fingerprint'  serves  as  an  indicator  of  epithelial  cancers 
(cancers  in  the  area  near  the  surface  of  various  tissues  such  as  skin,  breast,  prostate  etc).  This 
article  describes  Teraview's  (University  of  Leeds)  imaging  system,  its  development  and 
challenges  to  be  overcome. 

11.  O'Hara,  J.;  Grischkowsky,  D.;  “Coherent  THz  image  processing  and  enhancement” 
Conference  on  Lasers  and  Electro-Optics,  2003.  CLEO  '03.,  1-6  June  2003  Pages:3  pp. 

This  report  uses  image  processing  to  achieve  spatial  resolution  enhancement  for  a  coherent, 
synthetic,  phased-array  THz  imaging  system.  The  article  states  that  these  techniques  eliminate 
higher  order  scattering  processes,  or  clutter,  and  recover  otherwise  lost  resolution  in  images  of 
two  abutting  397  /spl  mu/m  diameter  spheres. 

12.  Zhang,  X.-C.;  “Recent  progress  of  terahertz  imaging  technology”  Conference  on 
Optoelectronic  and  Microelectronic  Materials  and  Devices,  1 1-13,  Pages:  1-6,  Dec.  2002. 

A  perspective  on  state  of  art. 

13.  Wallace,  V.P.;  Arnone,  D.A.;  Woodward,  R.M.;  Pye  RJ;  “Biomedical  applications  of 
terahertz  pulse  imaging”,  [Engineering  in  Medicine  and  Biology,  2002.  24th  Annual 
Conference  and  the  Annual  Fall  Meeting  of  the  Biomedical  Engineering  Society] 
EMBS/BMES  Conference,  2002.  Proceedings  of  the  Second  Joint  ,  Volume:  3  ,  23-26 
Pages:2333  -2334,  Oct.  2002 

Explains  the  biomedical  opportunities  for  THz  imaging  since  the  technology  behind  optically- 
driven  terahertz  (THz)  system  is  now  well  developed  and  commercial  applications  of  this  new 
technology  are  now  beginning  to  emerge.  The  study  describes  the  application  of  terahertz 
pulse  imaging  for  imaging  diseased  tissue  over  a  region  of  the  electromagnetic  spectrum,  the 
THz  region  which  has  only  recently  opened  up  for  this  kind  of  investigation. 
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14. 


Wallace,  V.P.;  Cole,  B.C.;  Woodward,  R.M.;  Pye,  R.J.;  Arnone,  D.A.;  “Biomedical 
applications  of  terahertz  technology”  Lasers  and  Electro-Optics  Society,  2002.  LEOS  2002. 
The  15th  Annual  Meeting  of  the  IEEE  ,  Volume:  1  ,  Pages:308  -  309,  10-14  Nov.  2002  . 
Similar  to  13. 

15.  Herrmann,  M.;  Tani,  M.;  Watanabe,  M.;  Sakai,  K.;  “Terahertz  imaging  of  objects  in 
powders”  Optoelectronics,  IEE  Proceedings- ,  Volume:  149  ,  Issue:  3  ,  Pages:  116  -  120,  June 
2002,. 

The  use  of  terahertz  radiation  for-detecting  and  imaging  objects  hidden  in  powders  is 
demonstrated  in  this  paper.  A  simple  method  for  distinguishing  between  objects  of  different 
materials  even  if  their  thickness  is  unknown  is  explained  in  this  report.  The  potential  of  THz 
radiation  for  measuring  the  local  density  profile  of  powders  is  also  demonstrated. 

16.  Ferguson,  B.;  Wang,  S.H.;  Abbott,  D.;  Zhang,  X.-C.;  “T-ray  tomographic  imaging”  IEEE 
Tenth  International  Conference  on  Terahertz  Electronics  Proceedings,  2002.,  9-10  Sept.  2002. 
Reports  about  three  THz  wave  tomographic  imaging  modalities:  T-ray  Computer 
Tomography  (T-ray  CT),  T-ray  Diffractive  Tomography  (T-ray  DT)  and  tomographic 
imaging  with  Fresnel  lens.  T-ray  CT  provides  3D  images  and  reveals  information  about  the 
target  for  functional  imaging.  T-ray  DT  achieves  sub-mm  spatial  resolution  and  binary  lens 
leads  to  novel  imaging  modalities  with  broad  band  T  ray  pulses. 

17.  Zhang,  X.-C.;  “Tomographic  imaging  with  terahertz  pulses”  Lasers  and  Electro-Optics 
Society,  2002.  LEOS  2002.  The  15th  Annual  Meeting  of  the  IEEE  Volume:  2,  Pages:  528  - 
529,  10-14  Nov.  2002. 

Explains  the  advantages  and  the  procedure  involved  in  using  Fresnel  lens  for  tomographic 
imaging  with  THz  radiation.  The  linearly  frequency-dependent  focal  length  of  a  Fresnel  lens 
allows  tomographic  imaging  of  a  target  using  multiple  frequencies.  Objects  at  various 
locations  along  the  beam  propagation  path  are  uniquely  imaged  on  the  same  imaging  plane 
using  a  Fresnel  lens  with  different  frequencies  of  the  imaging  beam.  This  procedure  allows 
the  reconstruction  of  an  object's  tomographic  contrast  image  by  assembling  the  frequency- 
dependent  images,  and  provides  a  new  tomographic  imaging  modality. 

18.  Chamberlain,  J.M.;  Zinovev,  N.N.;  Fitzgerald,  A.J.;  Berry,  E.;  Walker,  G.C.;  Handley,  J.W.; 
Smith,  M.A.;  Miles,  R.E.;  “What  constitutes  a  useful  and  meaningful  terahertz  image?” 
IEEE  Tenth  International  Conference  on  Terahertz  Electronics  Proceedings,  2002.  Pages:  81  - 
84,  9-10  Sept.  2002. 
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Attempts  to  answer  the  question  how  are  the  THz  images  best  processed,  analyzed,  presented, 
classified  ad  interpreted  with  a  clinical  setting?  Computer  vision  and  data  compression 
techniques  are  also  discussed  in  the  process.  It  concludes  that  compounded  images,  involving 
pulse  broadening  as  an  image  parameter,  are  highly  promising  routes  to  visualization. 

19.  Zinov'ev,  N.N.;  Fitzgerald,  A.F.;  Strafford,  S.M.;  Wood,  D.J.;  Carmichael,  F.A.;  Miles,  R.E.; 
Smith,  M.A.;  Chamberlain,  J.M.;  “Identification  of  tooth  decay  using  terahertz  imaging 
and  spectroscopy”  Infrared  and  Millimeter  Waves,  2002.  Conference  Digest.  Twenty 
Seventh  International  Conference  on  ,  22-26  Sept.  2002  Pages:13  -  14  (September  2002). 

The  terahertz  (THz)  frequency  spectroscopic  imaging  studies  of  teeth  are  reported.  The  aim  of 
the  experiment  was  to  establish  the  characteristic  properties  of  the  enamel  and  dentine  at  these 
high  frequencies.  Changes  to  the  THz  characteristics  as  a  result  of  various  types  of  tooth 
decay  are  reported  showing  the  potential  of  this  technique  for  dental  diagnosis. 

20.  Herrmann,  M.;  Tani,  M.;  Sakai,  K.;  Watanabe,  M.;  “Towards  multi-channel  time-domain 
terahertz  imaging  with  photoconductive  antennas”  International  Topical  Meeting  on 
Microwave  Photonics,  2002.,  Pages:317  -  320,  5-8  Nov.  2002 

Explains  the  working  and  results  of  a  small  and  cheap,  yet  efficient  electronic  circuit  for 
reading  out  THz  signals  from  photoconductive  antennas.  An  8-channel  version  of  this  circuit 
has  been  used  for  THz  imaging  by  simultaneously  recording  sets  of  8  time-domain  THz 
waveforms  with  an  array  of  8  photoconductive  antennas. 

21.  Loffler,  T.;  Bauer,  T.;  Siebert,  K.J.;  Roskos,  H.G.;  Fitzgerald,  A.;  Czasch,  S.;  “Investigation 
of  tumor  recognition  by  Terahertz  dark-field  imaging”  Lasers  and  Summaries  of  Papers 
Presented  at  the  Electro-Optics,  2002.  CLEO  '02.  Technical  Digest.  19-24  May  2002  . 
Demonstrates  that  the  THz  dark-field  imaging,  which  allows  enhancing  of  imaging  contrast 
especially  at  tissue  boundaries,  holds  for  the  clinical  distinction  between  beginning  and 
malignant  tumors  as  this  is  usually  based  on  the  differences  in  the  structure  of  the  tumor 
boundary.  This  article  states  that  dark-field  imaging  deals  with  contrast  enhancement  by  the 
detection  of  that  part  of  the  radiation  which  is  deflected  out  of  the  beam-propagation  direction 
by  either  diffraction  or  scattering  in  the  sample,  and  requires  blocking  of  the  ballistic  part  of 
the  radiation. 

22.  Mittleman,  D  M.;  “T-ray  imaging:  new  possibilities  in  the  far  infrared”  Microwave 
Symposium  Digest,  2002  IEEE  MTT-S  International  ,  Volume:  3  ,  Pages:1575  -  1578,  2-7 
June  2002 
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Describes  one  technique  in  terahertz  imaging,  which  permits  time-of-flight  imaging  with 
depth  resolution  well  below  the  Rayleigh  limit.  These  include  implementations  in  reflection 
geometry,  where  the  short  duration  of  the  THz  pidses  permits  three-dimensional  image 
reconstruction  using  time-of-flight  techniques.  The  paper  also  discusses  a  new  imaging 
technique  involving  interferometry,  which  can  be  used  to  resolve  two  reflecting  surfaces 
which  are  separated  by  much  less  than  the  coherence  length  of  the  radiation. 

23.  Foulds,  A.P.;  Chamberlain,  J.M.;  Berry,  E.;  “Terahertz  imaging  of  wound  healing”  High 
Frequency  Postgraduate  Student  Colloquium,  2002.7th  IEEE  ,  8-9  Sept.  2002  . 

Explains  an  experimental  system  which  was  being  built  to  produce  images  in  the  little  used 
frequency  range  of  0.1-3  THz  of  the  electromagnetic  system.  A  system  is  described,  based  on 
a  femtosecond  laser  using  optical  rectification  to  generate  and  electro-optic  sampling  to  detect 
the  terahertz  radiation.  Different  types  of  THz  images  are  introduced,  and  the  application  of 
THz  imaging  to  wound  healing  is  discussed. 

24.  Schildknecht,  C.;  Kleine-Ostmann,  T.;  Knobloch,  P.;  Rehberg,  E.;  Koch,  M.;  “Numerical 
image  enhancement  for  THz  time-domain  spectroscopy”  IEEE  Tenth  International 
Conference  on  Terahertz  Electronics  Proceedings,  2002.,  Pages:  157  -  160,  9-10  Sept.  2002. 
Explains  the  application  of  numerical  image  enhancement  methods  to  THz  time-domain 
spectroscopy  images  in  order  to  improve  the  quality  of  experimental  data  suffering  from  low 
resolution  and  finite  focus  diameter  aperture  effects.  As  an  example,  the  level  of  detail  in  a 
biomedical  image  of  a  histopathological  sample  is  shown  to  increase  significantly.  The  article 
presents  the  Jansson-Van-Cittert  algorithm  as  a  possible  approach  to  numerical  image 
enhancement  of  THz  time  domain  spectroscopic  images.  Its  states  that  numerical  image 
processing  has  more  applications  in  THz  time  domain  spectroscopy  than  image  enhancement. 
It  also  states  that  digital  filtering  can  be  used  to  calculate  the  derivative  of  the  image  profile  in 
order  to  reveal  and  recognize  edges  enabling  applications  in  THz  packaging  inspection. 

25.  Schmidt,  L.-P.;  Biber,  S.;  Rehm,  G.;  Huber,  K.;  “THz  measurement  technologies  and 
applications”  14th  International  Conference  on  Microwaves,  Radar  and  Wireless 
Communications,  2002.  MIKON-2002.,  Volume:  2  ,  Pages:581  -  587,  20-22  May  2002. 

Looks  at  THz  detectors. 

26.  Herrmann,  M.;  Tani,  M.;  Sakai,  K.;  Watanabe,  M.;  “Multi-channel  signal  recording  with 
photoconductive  antennas  for  THz  imaging”  IEEE  Tenth  International  Conference  on 
Terahertz  Electronics  Proceedings,  2002.,  9-10  Sept.  2002  Pages:28  -  31 
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Similar  to  [20].  Explains  the  working  and  results  of  a  small  and  cheap,  yet  efficient  electronic 
circuit  for  reading  out  THz  signals  from  photoconductive  antennas.  An  8-channel  version  of 
this  circuit  has  been  used  for  THz  imaging  by  simultaneously  recording  sets  of  8  time-domain 
THz  waveforms  with  an  array  of  8  photoconductive  antennas. 

27.  Shaohong  Wang;  Ferguson,  B.;  Mannella,  C.;  Abbott,  D.;  Zhang,  X.-C.;  “Powder  detection 
using  THz  imaging”  Summaries  of  Papers  Presented  at  the  Quantum  Electronics  and  Laser 
Science  Conference,  2002.  QELS  '02.  Technical  Digest.,  19-24  May  2002 

Contains  a  brief  description  of  the  possibilities  of  using  terahertz  waves  in  detecting  harmful 
substances  such  as  explosives,  illicit  drugs  and  bioterrorism  substances  such  as  the  Anthrax 
causing  bacterium  Bacillus  anthracic.  This  article  states  that  the  nondestructive  scanning  of 
mail  envelopes  and  packages  by  a  chirped  probe  beam  THz  imaging  system  is  a  topic  of 
current  research  interest. 

28.  Kleine-Ostmann,  T.;  Knobloch,  P.;  Koch,  M.;  Hoffmann,  S.;  Hofmann,  M.;  Hein,  G.;  Pierz, 
K.;  “Compact  and  cost-effective  continuous  wave  THz  imaging  system”  Summaries  of 
Papers  Presented  at  the  Lasers  and  Electro-Optics,  2002.  CLEO  '02.  Technical  Digest., 
Pages:405  -  406,  19-24  May  2002. 

THz  radiation  is  generated  via  photo  mixing  in  a  LT-GaAs  dipole  antenna.  Research  shows 
that  during  the  last  years  THz  technology  has  been  developed  to  a  stage  where  first 
commercial  systems  are  available.  Yet,  the  THz  imaging  systems  developed  so  far  all  work  in 
the  time-domain  using  THz  pulses.  Hence,  they  rely  on  expensive  and  bulky  femtosecond 
lasers  as  a  key  component.  Therefore  this  article  details  a  terahertz  imaging  setup  that  can 
generate  frequencies  from  0.2  to  0.5  THz. 

29.  Siebert,  K.J.;  Quast,  H.;  Leonhardt,  R.;  Loffler,  T.;  Thomson,  M.;  Bauer,  T.;  Roskos,  H.G.; 
Czasch,  S.;  “All-optoelectronic  CW  THz-imaging”  Summaries  of  Papers  Presented  at  the 
Lasers  and  Electro-Optics,  2002.  CLEO  '02.  Technical  Digest.,  Pages:635  -  636,  19-24  May 
2002  . 

Demonstrates  CW  THz  imaging  with  complex  sample  and  compares  the  results  with  images 
taken  with  a  state-of-the-art  pulsed  system  based  on  an  Ti:sapphire  amplifier  laser  operating  at 
1-kHz  repetition  rate. 

30.  T.  D.  Dorney,  W.  W.  Symes,  R.  G.  Baraniuk,  and  D.  M.  Mittleman,  “Terahertz  multistatic 
reflection  imaging,”  Journal  of  the  Optical  Society  of  America  A,  19,  2002.. 

Applies  seismic  migration  type  techniques  to  image  conducting  objects  in  dielectrics. 
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31. 


Shen,  Y.C.;  Upadhya,  P.;  Davies,  A.G.;  Linfield,  E.H.;  “Light-induced  difference 
Terahertz  spectroscopy  and  its  biomedical  applications”  Terahertz  Electronics 
Proceedings,  2002.  IEEE  Tenth  International  Conference  on  ,  Pages:  161  -  164,  9-10  Sept. 
2002. 

Uses  laser- induced  difference  TEIz  spectroscopy  to  investigate  three  samples  with  different 
lifetimes.  The  spectroscopy  system  is  based  on  a  10  nJ  titanium  sapphire  laser  with  a  pulse 
duration  of  12  fs  and  a  centre  wavelength  of  790  nm.  A  lifetime  of  about  50  ps  has  been 
determined  for  the  SI-GaAs  sample,  whereas  the  lifetime  of  the  EIR-Si  samples  was  found  to 
be  much  larger  than  the  time  interval  between  two  successive  laser  pulses  (12  ns).  As  a  result, 
the  differential  THz  signal  is  about  twenty  times  larger  than  that  for  SI-GaAs.  This  article  also 
observed  that  the  TEIz  pulse  arrives  at  the  detector  100  fs  earlier  when  it  transmitted  through 
an  optically  excited  EIR-Si  wafer. 

32.  De  Lucia,  F.C.;  “THz  spectroscopy  -  techniques  and  applications”  Microwave  Symposium 
Digest,  2002  IEEE  MTT-S  International,  Volume:  3,  Pages:  1579  -  1582,  2-7  June  2002. 
Reviews  established  techniques  and  discuss  new  approaches  in  THz  spectroscopy. 
Spectroscopy  and  elated  applications  have  accounted  for  a  large  proportion  of  the  work  in  the 
THz  spectral  region.  These  have  included  laboratory  studies  of  the  radiative  and  collisional 
properties  of  molecules,  ions,  and  radicals,  as  well  as  field  applications  such  as  atmospheric 
remote  sensing  and  molecular  radio  astronomy. 

33.  Kadlec,  R.;  Nemec,  H.;  Kuzel,  P.;  “Frequency-domain  approach  to  evaluation  of  data 
obtained  in  optical  pump-terahertz  probe  experiments”  Twenty  Seventh  International 
Conference  on  Infrared  and  Millimeter  Waves,  2002.  Conference  Digest.  ,  Pages:23  -  24  22- 
26  Sept.  2002. 

Describes  the  increasing  effort  to  develop  optical  pump-terahertz  probe  (OPTP)  experiments 
as  a  tool  for  studying  transient  processes  following  excitation  by  an  optical  pulse,  e.g.  in 
semiconductors,  superconductors.  A  method  is  necessary  to  find  the  unknown  non¬ 
equilibrium  part  of  the  response  function  so  the  authors  of  this  article  have  developed  a 
general  frequency  domain  theory  which  can  serve  as  a  guide  for  performing  the  experiment. 

34.  Chen,  Q.,  Zhang,  X.-C.,  "Semiconductor  dynamic  aperture  for  near-field  terahertz  wave 
imaging",  Selected  Topics  in  Quantum  Electronics,  IEEE  Journal  of  Quantum  Electronics, 
Volume:  7  ,  Issue:  4  ,  pages  608  -  614,  July-Sept.  2001. 

Describes  and  discusses  near-field  terahertz  wave  imaging  with  a  dynamic  aperture  created  on 
a  semiconductor  wafer.  The  spatial  resolution  of  a  near-field  terahertz  wave  imaging  system 
with  a  dynamic  aperture  created  on  a  GaAs  wafer  is  determined  by  the  diameter  of  the  gating- 


40 


beam- induced  thin  photo  carrier  layer.  With  a  dynamic  aperture  created  on  a  GaAs  wafer,  the 
experiment  achieved  a  sub  wavelength  spatial  resolution  of  -A/10.  In  addition,  near-field 
terahertz  wave  imaging  with  a  dynamic  aperture  overcomes  the  limitations  in  the  reported 
methods  of  near-field  terahertz  wave  imaging.  Various  terahertz  images  are  presented  to  show 
the  possible  applications. 

35.  Johnson,  J.L.;  Dorney,  T.D.;  Mittleman,  D.M.;  “Interferometric  imaging  with  terahertz 
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